The probability distribution of the conductance in anisotropic systems 
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We investigate the probability distribution p{g) of the con- 
ductance g in anisotropic two-dimensional systems. The scal- 
ing procedure applicable to mapping the conductance distri- 
butions of localized anisotropic systems to the corresponding 
isotropic one can be extended to systems at the critical point 
of the metal-to-insulator transition. Instead of the squares 
used for isotropic systems, one should use rectangles for the 
anisotropic ones. At the critical point, the ratio of the side 
lengths must be equal to the squre root of the ratio of the criti- 
cal values of the quasi-one-dimensional scaling functions. For 
localized systems, the ratio of the side lengths must be equal 
to the ratio of the localizaticm lengths. 

The presence of disorderly may allow a system to make 
a transition from metallic to insulating behaviour by 
varying the Fermi energy in an energy range where both 
extended and localized states are found, separated by a 
mobility edge. Characterizing this transition, one can 
employ transport properties, such as the conductance, or 
properties of the system's eigenstates, such as the corre- 
lation length for extended, metallic states or the local- 
ization length ^ for insulating states. At the mobility 
edge, a determination of the complete probability distri- 
bution p{g) of the conductance g (in units of e'^/h) is 
needed. The critical point of the transition from metal- 
lic states to Anderson localized onesB is of particular 
interest. The distributions are well-known to be nor- 
mal and log-normal off the mobility edge towards the 
extended and the localized regime respectively, whereas 
the exact forHi|j3f the critical distribution is still under 
investigation.BiiEl For example, contrary to expectations, 
the critical distribution seems to vary even within the 
same universality class, depending on the boundary con- 
ditionsjperpendicular to the direction in which transport 
occursJjia Also, questions about the exact form of the 
large-g tail {g > 1) remain unanswered. Where calcula- 
tions in 2 -t- e dimensionsij indicate higher cumulants to 
diverge with system_size, leading to a power law tail, nu- 
merical calculationaO in three dimensioiis and analytical 
results for quasi-one-dimensional wiregj show an expo- 
nential decay. 

Anisotropic systerpa Jiave recently been the focus of 
particular attention.Eirt3 It is generally accepted that 
anisotropy does not change the universality class and that 
isotropic results can be recovered by performing a proper 
scaling of the anisotropic results. For anisotropic sys- 
tems in a localized state, it is reasonable to assume that 
scaling the dimensions of the system by the correspond- 
ing localization lengths will make the system effectivdx 
isotropic. This procedure has been applied successfuUytll 
to the scaling function A ~ Xm/M, which is a function of 



£,/M, where Am denotes the finite size localization length 
of a quasi-one-dimensional strip of finile width M and 
Am — > C 8.S M ^ oo. It was also shownEZi that the same 
scaling procedure works for the probability distribution 
in such a system. 

In order to test the approach for critical states one 
must either face the numerical challenge of large three- 
dimensional systems or take into account additional in- 
teractions (bejpiid the disorder potential) such as spin- 
orbit coupling.li3 Another possibility is the introduction 
of exterual magnetic fields as e.g. in integer quantum Hall 
systemstj|-ar tight binding models with random mag- 
netic fluxBj However, as a result of Anderson localiza- 
tion, extended states do not exist in two-dimensional 
systems of non-interacting electrons in a magnetic field, 
except at a singular energy near the center of each of 
the Landau subbands. At these critical energies Ec the 
localization length ^ diverges with a critical exponent v: 

Because significant finite size effects have to be ex- 
pected, we decided to concentrate our research on two- 
dimensional systems although the exact form of the crit- 
ical distribution of the conductance depends on the di- 
mensionality of the system.Q The investigation of the 
self-averaging quantity ^ in integer qjuantum Hall sys- 
tems yielded very encouraging results )l3 supporting the 
expectation that quantities of anisotropic systems can in- 
deed be mapped to isotropic values by a simple rescaling 
scheme. 

In this paper, we show a method of mapping the proba- 
bility distributions of the conductance of anisotropic two- 
dimensional systems with a magnetic field perpendicular 
to the plane or with spin-orbit coupling to the probabil- 
ity distribution of the conductance for the corresponding 
isotropic system at the critical point, using a tight bind- 
ing model. It turns out that the ratio of the squares of 
the side lengths L^, Ly of the anisotropic system should 
be chosen equal to the ratio of the critical values AJ^ 
of the quasi-one-dimensional scaling functions: 
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In the following, we first describe the models and the 
numerical method we employed. Then we present and 
discuss our numerical results and finally summarize the 
conclusions of this work. 

The tight binding model uses the Hamiltonian 
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where n, n' denotes the lattice site. Without spin-orbit 
interaction the "variables" r, r' take on only one value 
and the hopping integrals Vn,n' are scalar, otherwise they 
are 2x2 matrices and the spin variables take on the 
values 1 or —1. In either case the site energies £„ are 
independent of t and we take into account interactions 
only between neighbouring lattice sites. 

An external magnetic field enters the Hamiltonian via 
its vector potential A (V x A = B), which appears in 
the phases of the hopping integrals: 



_ -27rj(e//i) f^^' A(r)dr 
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The integral connects the lattice sites n (at r„) and n' 
(at r„/) in a straight line. For the sytems under consider- 
ation, where the magnetic induction B is perpendicular 
to the plane of the two-dimensional lattice, the gauge for 
the vector potential can be chosen such that the phases 
vanish in the direction perpendicular to A and are integer 
multiples of some number 27ra in the direction parallel 
A. The value of the parameter a then completely charac- 
terizes the influences of the magnetic field on the system. 
For rational a, the denominator determines the number 
of bands in the density of states of the system without 
disorder. p. 

The Evangelou-Ziman modeled incorporates spin- 
orbit coupling by using the following hopping integrals: 



6r,r' + fJ.iJ2K,T'tn, 
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where v = x,y,z and a'^ are the Pauli matrices. The 
parameter fi characterizes the strength of the spin-orbit 
interaction. 

Both systems may be made anisotropic by chosing the 
value of to be different in the two directions within the 
plane. Otherwise this parameter is a constant, indepen- 
dent of lattice site n. We bring disorder to the system by 
chosing all the site energies independently from a rect- 
angular distribution of width W centered at 0, so that 
is a measure of the strength of the disorder. The pa- 
rameters are also randomly selected from a uniform 
distribution on [—1/2,1/2]. The energy scale is set by 
the larger of the two values for which is taken to be 
unity. 

We c^culate the conductance from the Landauer 
formulaca 
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where t is the transmission matrix. We suppose two 
semi-infinte leads are attached to opposite sides of the 
sample. Then t determines the transmission of an 
electron through the sample. The numerical proce- 
dure is |based on the algorithms published by Ando and 
Pendry£3 

The critical conductance distributions we calculated 
for isotropic quantum Hall systems at different disorder 



strengths show that finite size effects become stronger 
the weaker the disorder. Where systems with W — 4.0 
and W = 2.0 show a basically size-independent critical 
distribution of the conductance for squares of 64 x 64 
lattice sites, at = 0.5 finite size effects are still some- 
what noticeable up to systems with 192 x 192 lattice 
sites. The anisotropic quantum Hall systems we inves- 
tigated are characterized by a = |. The anisotropics 
we chose were t^/t^ = 0.5 at = 0.5 and t°/t° = 0.8 
at W = 0.1. The latter was chosen mainly because we 
already had the data for the quasi-one-dimensional scal- 
ing function. As the disorder is even weaker than in the 
first case, finite size effects are even stronger, and even at 
240 X 240 lattice sites the conductance distribution is far 
from the one we expect from our calculations of isotropic 
systems. Therefore, we will not be able to show that our 
procedure maps the two anisotropic conductance distri- 
butions to the critical distribution of isotropic systems 
for this extreme case. We will however be able to prove 
the somewhat weaker claim that our method transforms 
the two anisotropic distributions so that both have the 
same shape. In a square system, one expects that the 
distribution in the difficult hopping direction shows a 
more localized character than the one in the easy hop- 
ping direction. In an isotropic system, the distribution 
obviously cannot depend on the direction of transport. 
By making the system rectangular rather than square, 
i.e. shorter in the difficult hopping direction, it should 
be possible to obtain distributions in the two directions 
that are the same, thus making the anisotropic system 
effectively behave isotropically. 

The task now is "How do we choose the correct ra- 
tio of side lengths ofj-Lha, rectangle?" From the research 
on localized systemsllifEZi we know that in those cases, 
the ratio should be equal to the ratio of the localization 
lengths: 
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for localized systems, as these are obviously the appro- 
priate length scales in their respective directions. This is 
of no use for critical systems as both localization lengths 
diverge at the transition. A closely corresponding non- 
diverging quantity is, however, available in the scaling 
function Am = \m/M, which has a finite critical value, 
independent of the system width M . The finite size local- 
ization lengths Xm,x and \M,y have and respectively 
as their large-M limits, and for big enough systems we 
can approximate Eq. y by 
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for localized systems at "large enough" M . The meaning 
of M in this context would be that of the system width 
perpendicular to the direction in which the localization 
length is measured, i.e. M — Ly in Xm,x and M — in 
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^M,y Now, by multiplying both sides of Eq. |^ by L^/Ly 
we have 
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for "large" localized systems. Now Am is a continu- 
ous function of E and for "large enough" systems at Ec 
should have reached its critical value. Therefore, we ar- 
rive at the conclusion that 
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should be the correct ratio for critical systems in order 
to make th|aBi,^have like isotropic ones. Noting the 
relationshipQilS'tZl 



A'^ • A'^ = A? 
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we can write Eq. ^ in the alternate form 
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We tested this prediction on the system with t^/t^ — 
0.5 and W = 0.5, where the ratio y^Aj/A^ is roughly 1.5. 

The result is shown in Fig. |l] together with the critical 
distribution for the isotropic system. The agreement is 
very good. For the other system with = 0.8 and 

W = 0.8, we have to deal with stronger finite size effects 
and cannot expect to approach the form of the critical 
distribution we see in Figs. |l| for reasonable system sizes. 
Instead we merely show in Fig. |^ how the critical distri- 
butions change with the ratio of side lengths. The best 
value for the ratio according to Eq. ^ would be roughly 
1.23. Fig. ^ shows results for ratios of 1.0, 1.25 and 1.5. 
The averages of \n{g) for the easy hopping direction de- 
crease with increasing ratio from —4.09 for the square 
to —4.40 and —5.30, while the averages for the difficult 
hopping direction increase from —4.94 for the square to 
—4.18 and —4.03. Similarly, the standard deviations in- 
crease for the easy hopping direction from 1.90 for the 
square to 2.07 and 2.30, while they decrease for the dif- 
ficult hopping direction from 2.24 for the square to 2.01 
and 1.91. The values for a ratio of 1.25 are not equal 
but reasonably close, so that for larger systems, where 
a ratio of 1.23 might be practicable, we expect a better 
agreement of the two probability distributions. 

Taking the best-ratio rectangle as the "undeformed" 
base, we can also see from Fig. ^ that similar "deforma- 
tions" have similar effects in the two directions, that is, 
reducing the ratio by a factor 7 in one direction will cause 
the ensemble average < \D.{g) > in that direction to in- 
crease, and the standard deviation to decrease, while the 
trend is opposite in the perpendicular direction. How- 
ever, by reducing the ratio by the same factor 7 in the 



other direction, will result roughly in the same distri- 
butions as before, but with the one associated with the 
easy direction before now assigned to the difficult direc- 
tion and vice versa. 

That the same procedure also works for sytems with 
spin-orbit coupling is shown in Fig. |^, where we plot 
the conductance distribution for an isotropic system with 
^ = 1.0 at i? = 0.1 and Wc = 6.7 together with that of 
two anisotropic systems, one with ^ — 1.0 and t'^/ty — 
0.1 at = 0.1 and Wc = 1.6, the other with ^ = 1.0 
and t°/t° = 0.2 a,t E ^ 0.1 and Wc = 2.6. The ratio 
of sidelengths, according to Eq. |9[ should be 23.0 for the 
latter. We chose 40 x 40 lattice sites for the isotropic 
system, 10 x 230 lattice sites for the strongly anisotropic 
one and 20 x 185 lattice sites for the anisotropic system 
with the weaker anisotropy. Again the agreement is very 
good. 

We have shown that the scaling procedure applica- 
ble to mapping the conductance distributions of local- 
ized anisotropic systems to the corresponding isotropic 
one can be extended in a straightforward manner to sys- 
tems at the critical point of the Anderson localization- 
delocalization transition in both unitary and symplectic 
two-dimensional systems. Instead of the squares used 
for isotropic systems, one should use rectangles for the 
anisotropic ones, with a ratio of side lengths equal to 
the square root of the ratio of the critical values of the 
quasi-one-dimensional scaling function. 
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FIG. 1. The conductance distributions of an anisotropic 
rectangular system with a ratio of side lengths chosen accord- 
ing to Eq. 1^. For comparison the corresponding distribution 
of an isotropic system is shown as well. 
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FIG. 2. Conductance distributions of an anisotropic sys- 
tem for varying ratios of the side lengths. The left panels refer 
to transport in the easy hopping direction, the right panels to 
transport in the difficult hopping direction. The ratio of side 
lenghts is 1.0 for the top row, 1.25 for the middle row and 1.5 
for the bottom row. 
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FIG. 3. The conductance distributions for an isotropic 
and two anisotropic systems with spin-orbit coupling. 
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